Strongly interacting Fermi systems in expansion: 
Prom cold atoms to color superconductivity 
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We investigate the expansion proposed recently as a strategy to include quantum fluctuation 
effects in the nonrelativistic, attractive Fermi gas at and near unitarity. We extend the previous 
results by calculating the next-to-leading order corrections to the critical temperature along the 
whole BCS-BEC crossover. We demonstrate explicitly that the extrapolation from the mean-field 
approximation, based on the expansion, provides a useful approximation scheme only on the 
BCS side of the crossover. We then apply the technique to the study of strongly interacting rela- 
tivistic many-fermion systems. Having in mind the application to color superconductivity in cold 
dense quark matter, we develop, within a simple model, a formalism suitable to compare the effects 
of order parameter fluctuations in phases with different pairing patterns. Our main conclusion is 
that the relative correction to the critical temperature is to a good accuracy proportional to the 
mean-field ratio of the critical temperature and the chemical potential. As a consequence, it is 
significant even rather deep in the BCS regime, where phenomenologically interesting values of the 
quark-quark coupling are expected. Possible impact on the phase diagram of color-superconducting 
quark matter is discussed. 
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I. INTRODUCTION 

Strongly interacting many-fermion systems have been 
a theoretical challenge for decades. While the Bardeen- 
Cooper-Schrieffer (BCS), or mean- field (MF), theory 
provides an excellent description of conventional metallic 
superconductors, it is still inappropriate for a large range 
of other systems, from high-temperature superconductors 
to nuclear (and more recently, quark) matter. 

Great theoretical interest was triggered by the con- 
jecture that when the interaction strength is increased, 
BCS-type superconductivity evolves smoothly to the 
Bose-Einstein condensation (BEC) of tightly bound 
difermion molecules P, 0, Q. The crucial observation 
in this respect was that the standard BCS superconduct- 
ing ground state has the same form as the ground state 
of a condensed Bose gas once the composite operator cre- 
ating a Cooper pair is identified with that of a bosonic 
quasiparticle. However, in spite of the successful unified 
description of BCS superconductivity and BEC, a quan- 
titative understanding of the crossover between the two 
regimes was missing. 

Spectacular progress in this direction has been made in 
the past decade thanks to the experiments using ultracold 
atomic Fermi gases [1, 0]- From the theorists' point 
of view, these provide an ideal tool to test the developing 
many-body techniques. In particular, a lot of interest has 



been attracted by the transition regime between the BCS 
and BEC limits — the unitary Fermi gas. In this case, the 
two-body scattering length is much larger (ideally infi- 
nite) than any other characteristic scale of the system 
such as the interaction range and interparticle distance. 
This, on the one hand, leads to intriguing universal be- 
havior with connections to other branches of physics such 
as the quark-gluon plasma or even string theory . 
On the other hand, it poses a challenging problem due 
to the lack of a small expansion parameter. 

Several approaches have been suggested to deal with 
the Fermi gas at large scattering length, including vari- 
ous self-consistent resummation techniques for the many- 
body Green's functions [l^, |Tl[ or the scattering ma- 
trix [l2l|. e xpansion in the dimensionality of the space 
[T^. Il4l. Iisj. or the 1/iV expansion [H, An early re- 
view of the many-body approaches to the crossover prob- 
lem may be found in Ref. [l^. The predictions of the 
analytic approximation schemes as well as the available 
experimental results are now bein g te sted by increasingly 
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precise numerical simulations 

Also in high-energy physics has the mechanism of 
Cooper pairing proven extremely fruitful. The Nambu- 
Jona-Lasinio (NJL) model [l^, , constructed in direct 
analogy with the BCS theory of superconductivity, was 
one of the first models of dynamical symmetry break- 
ing. For its simplicity and generality, it still remains a 
popular low-energy effective description of strongly inter- 
acting quark and nuclear matter, this article not being 
an exception. 

While the original NJL model dealt with dynami- 
cal breaking of chiral symmetry by particle-antiparticle 
correlations, the true Cooper pairing of two relativistic 
fermions near their Fermi surface appears in dense quark 
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matter, leading: to the so-called color superconductivity 
[25I [26| (see [l^] for a recent review). Due to the strong- 
interaction nature of quantum chromodynamics (QCD), 
quark matter at moderate densities is a typical example 
of a relativistic many-fermion system where a departure 
from the BCS-like behavior is to be expected. The early 
work in this respect focused on the structural change of 
Cooper pairs at strong coupling and precursor phenom- 
ena above the critical temperature for the superconduct- 
ing phase transition, in particular the appearance of the 
pseudogap in the spectrum [H, [H, [s^, [U [H, HI] . 

In principle, the size of the QCD running coupling de- 
pends just on the energy/momentum scale and is thus 
fixed by the density of the quark matter. However, 
in order to better understand the strong-coupling ef- 
fects in color superconductors, one often considers mod- 
els with variable coupling strength where the full BCS- 
BEC crossover can be studied O, [H, [H, [13, [H, [H, E^l . 
In the color superconductivity context, most of the cal- 
culations are still done using the MF approximation. 
The first attempts to include the fluctuations of the or- 
der parameter have adapted the Nozieres-Schmitt-Rink 
(NSRUheory [1, Hi, HI and the pseudo gap approxima- 
tion |4l| . commonly used in condensed-matter phy sics, 
and the Cornwall-Jackiw-Tomboulis formalism p3 . |43| , 
well known in high-energy physics. 

The goal of this paper is to investigate the applica- 
bility of the expansion to the study of strongly in- 
teracting Fermi gases. In Sec. |IT] we review the 
expansion for the nonrelativistic Fermi gas developed in 
Refs. [H, • We study the evolution of the next-to- 
leading-order corrections from the unitarity towards the 
BCS and BEC regimes and critically examine the virtues 
as well as shortcomings of the method. In Sec. IIIII we 
then apply the technique to relativistic superconductors, 
using a class of NJL-type models. In particular we esti- 
mate the correction to the critical temperature in color- 
superconducting quark matter and show that it is signifi- 
cant even for realistic values of the coupling strength. On 
account of the fact that the fluctuation effects are differ- 
ent for various competing phases, we propose a modifica- 
tion of the QCD phase diagram. Finally, in Sec. IIVI we 
summarize and conclude. 



II. NONRELATIVISTIC ATTRACTIVE FERMI 
GAS 

We consider here the gas of two nonrelativistic fermion 
species ("flavors") which will be referred to as ip^ and 
'0X- For simplicity they will be assumed to have equal 
masses and chemical potentials (and thus also densities) . 
Nevertheless, the results of the present analysis may be 
straightforwardly generalized to the case of a density im- 
balance, or, Fermi surface mismatch At low density, 
i.e., low characteristic momentum set by the Fermi scale, 
the two-body interaction is completely determined by the 
s-wave scattering length, a. 



Formalism 



In order to employ the expansion, one has to 

generalize the system by including N copies of the two 
fermion flavors. Following Refs. [la. [T7|. we write down 
the Euclidean Lagrangian in the form 
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N 
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N 

The sums over repeated indices will from now on be im- 
plicitly assumed. For iV = 1 the Lagrangian reduces to 
one describing a two-flavor gas with local contact attrac- 
tive interaction with strength g [43 |. In this extended 
version the coupling is rescaled as 5 — > g/N in order to 
ensure that the action scales naturally with N. 

Note that the Lagrangian ([T|) possesses, apart from the 
phase invariance generated by the total particle number, 
a symplectic symmetry, Sp(2iV). However, as will be- 
come clear soon, this symmetry remains unbroken by the 
Cooper pairing so that it does not give rise to any un- 
wanted Nambu-Goldstone (NG) bosons and hence does 
not affect the low-temperature thermodynamics of the 
system. 

As a next step the theory is bosonized by introducing 
the auxiliary field, (j) ~ '^V'iiV'iT' ^^'^ performing the 
Hubbard-Stratonovich transformation. The result is the 
nonlocal effective action. 
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where Q is the Nambu-space fermion propagator in pres- 
ence of the pairing field cf), 



-dr 



2m 



Eq. ^ can be interpreted as a classical action that de- 
fines a theory of the scalar field </> with self-interactions 
determined by the expansion of the action in powers of 
(/). The crucial observation made in Refs. [IB, [l3| is that 
since the action is proportional to N , the expansion of 
the partition function, or the thermodynamic potential, 
in powers of is equivalent to the expansion in loops. 
Here and in the following, the term "loop" is used to 
refer to a bosonic loop, unless explicitly indicated oth- 
erwise. Note that all loops containing fermions of the 
original theory ^ are resummed into the action ([2]), i.e., 
are included at the tree level with respect to bosons. 

Obviously, the leading order (LO) of the expan- 
sion is equivalent to the saddle-point approximation to 
the functional integral, that is, the usual MF approx- 
imation. The next-to-leading order (NLO) then incor- 
porates one-loop corrections, or, the Gaussian fiuctua- 
tions around the saddle point The expansion 
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thus provides a systematic ordering of the corrections to 
the MF approximation. We should nevertheless keep in 
mind that at the end of the calculation, we have to set 
N = 1. The way this extrapolation is performed is to 
be understood as a part of the definition of the method, 
which distinguishes it from other approaches with for- 
mally equivalent thermodynamic potential fsl, fisj . 

In general one calculates, in a given approximation 
scheme, the thermodynamic potential f2 as a function of 
the anticipated vacuum expectation value A of the field 
(j), and of the chemical potential fi. Their actual values 
in thermodynamic equilibrium are then determined by a 
simultaneous solution of the gap and number equations, 



(3) 



where n is the total particle density, related to the Fermi 
momentum by the usual expression, n = kp/Sn'^. It 
is well known that already at one-loop level, an attempt 
to solve the equations self-consistently leads to unphys- 
ical results, in particular the violation of the Goldstone 
theorem Veillette et al. [17] suggested to avoid 

this problem by an expansion of the gap and chemical 
potential simultaneously with the expansion of the ther- 
modynamic potential that follows from the action ([2]), 

17 = ivr!("' + + + . . . , 
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Comparing terms of the same order in the gap and num- 
ber equations Q , one obtains explicit expressions for the 
higher-order corrections to the MF values A^^\ In 
particular at NLO we find 



/id) 

A(i) 



(4) 



It is essential that all derivatives of the thermody- 
namic potential here are evaluated using the MF val- 
ues a'"), /z'"). One thus avoids the problems with self- 
consistency; in particular the NG boson of the sponta- 
neously broken symmetry is exactly masslcss [16]. 

When we merely wish to determine the critical tem- 
perature, the gap in Eq. ^ is fixed to zero and we solve 
for the temperature and chemical potential instead. The 
identity ^ then naturally modifies to one for the correc- 
tions of the variables of interest. 
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dAn<^'^ 
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We will comment later on the ambiguity that arises at 
this point, stemming from the fact that we can choose 
to solve Eq. ([3]) for 1/Tc (or any other function of Tc) 



and accordingly get a different type of 1/iV expansion — 
all derivatives with respect to T in Eq. ^ would simply 
turn into ones with respect to 1/T. 

Eqs. ([U and ([5]) of course have to be supplemented 
with an expression for the thermodynamic potential. For 
the sake of this section, we will need just the explicit form 
of the standard MF part. 
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and 



physical scattering length by 



2:^,, '^k = Ck - M, 
The bare coupling g is related to the 



using the usual notation ek = 2^^, Ck 



d^k 1 
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For further details, we refer the reader to the original 
literature [l^ [13] where all necessary formulas are de- 
rived in detail. In Sec. IIIII we will develop a relativistic 
formalism from which the present case will follow as a 
particular nonrelativistic limit. 



B. Numerical results 

1. Critical temperature 

While differing in other specific directions of investiga- 
tion, both Refs. [l^, [l7| addressed the question of spe- 
cial interest, the calculation of the critical temperature 
at unitarity. Veillette et al. obtained the result 



^ 0.4964 - i# 



(8) 



in units of the Fermi energy, ep — fcp/2m, whereas 
Nikolic and Sachdev calculated the correction to the in- 
verse temperature and got 



= 2.014- 



5.317 
N ■ 



(9) 



Both results are formally equivalent to order 1/A'' in the 
expansion, yet they yield dramatically different numbers 
when evaluated at iV = 1. Indeed, Eq. ([5]) becomes even 
meaninglessly negative and can be merely used to make 
the qualitative conclusion that the fluctuations decrease 
the critical temperature significantly. 

On the other hand, Eq. ([9]) leads to the critical temper- 
ature Tc — 0.14eF in remarkable agreement with the re- 
sult Tc = 0.152(7)eF, obtained by numerical Monte Carlo 
simulations [21], (2^ . However, it should be stressed that 
there is no a priori criterion that would tell us which 
observable to choose for the evaluation of the NLO cor- 
rection. One may think that it is the Lagrange multiplier 
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FIG. 1: Critical temperature and chemical potential as a func- 
tion of the inverse scattering length. Dashed lines: MF ap- 
proximation (LO in l/N). Solid lines: NLO calculation based 
on expansion of l/T. Dotted line in the upper panel: NLO 
value based on expansion of T as in Eq. ([5}. Dash-dotted 
lines: Self-consistent calculation using NSR theory. Dotted 
lines in the lower panel: First- and second-order perturbative 
approximations to the chemical potential, see Eq. (|10|l . 



gap equation. On the other hand, it then solves Eq. ([3]) 
self-consistently j6d| . 

In the BCS regime the critical temperatures calculated 
within the two approaches are consistent with each other, 
at least qualitatively. Around unitarity, NSR theory pre- 
dicts a well-known tiny maximum. On the contrary, the 
1 /N expansion gives a monotonic dependence on the in- 
verse scattering length, with the value at unitarity agree- 
ing very well with Monte Carlo calculations 21, 22], as 
noted above. However, as also emphasized there, the 
1/iV expansion should not be really trusted in this re- 
gion without an additional physical insight. Finally, in 
the BEG limit NSR theory converges to the expected 
critical temperature of the free Bose gas, Tc = 0.218eF. 
1/iV expansion fails to reproduce this asymptotic behav- 
ior. Although it cannot be seen in Fig. [1] its critical 
temperature keeps growing towards the BEG limit, al- 
beit with a decreasing slope. Indeed, one cannot really 
expect to get a constant asymptotics from a perturbative 
expansion around the rapidly increasing MF value. 

In the lower panel of Fig. [T] we also display the re- 
sults for the chemical potential at the critical tempera- 
ture. Here the most interesting is the evolution in the far 
BGS region. While in the MF approximation the chem- 
ical potential approaches its asymptotic value equal to 
the Fermi energy with an exponentially decreasing tail, 
upon including fluctuations the convergence turns into 
a much slower, power-law one. This is a purely pertur- 
bative effect, as clearly demonstrated by a comparison 
with the standard perturbative expansion for the chem- 
ical potential in the dilute Fermi gas f46] (the first- and 
second-order values are plotted using dotted lines) 



(3 ~ l/T rather than the temperature itself that is the 
natural variable of the thermodynamic potential. Still a 
deeper physical argument is obviously needed to resolve 
this ambiguity. Here we just remark that on the tech- 
nical level, it clearly arises from the truncation of the 
series for different observables; one cannot expect 
the expansion to be reliable when the NLO term is 
larger than the LO one [soj . 

In order to further study the size of the 1/iV corrections 
and the sensitivity to the choice of observable to evaluate 
them, we calculated the critical temperature at NLO as 
a function of the inverse scattering length, see Fig. [TJ 

Obviously, the evaluation of the critical temperature 
based on Eq. ([5]) (dotted line in the upper panel) is only 
reasonable in the far BGS regime (right end of the plot) 
where it roughly coincides with the 1/T-based approxi- 
mation (solid line). The l/T-based value will therefore 
be used exclusively in the following and will be referred 
to simply as the NLO critical temperature. For compar- 
ison we also show the critical temperature and chemical 
potential calculated using the NSR theory (dash-dotted 
line). NSR theory takes the fluctuation effects into ac- 
count only in the total particle number, but not in the 



f i^^i. , 4(11 -2 log 2) ^,^ 

— = 1 + — /cpaH TT^ (Kpa) H . (10) 

ep 6TT 157r^ 

Even though this formula holds only in the normal phase 
and at zero temperature, the agreement is excellent. The 
reason, of course, is that the pairing effects are suppressed 
exponentially and finite-temperature effects are also sup- 
pressed by the exponentially small value of the critical 
temperature in the BGS limit, so that they are both com- 
pletely negligible with respect to the perturbative correc- 
tions in Eq. (dni) m. 

The fact that the NLO in 1/N reproduces the pertur- 
bative expansion of the chemical potential up to second 
order can also be verified directly on the Feynman graph 
level. In Fig. [2] we show the perturbative contributions 
to the thermodynamic potential. In general the 1/N ex- 
pansion does not coincide with perturbation series — the 
diagram (d) is of second order in the interaction, yet only 
appears at the next-to-next-to-leading order in 1/7V. For- 
tunately, it turns out to vanish at zero temperature so 
that the 1/N expansion to NLO indeed contains the full 
second-order perturbative correction ^47j . 
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FIG. 2: Lowest-order perturbative contributions to the ther- 
modynamic potential of the Fermi gas. The free-gas graph (a) 
coincides with LO in The first-order graph (b) and the 

second-order graph (c) appear at NLO in \/N. The second- 
order graph (d) contributes only at NNLO. The dashed line 
denotes the bare propagator of the pairing field </>, g/N . 



2. BCS limit 

Since we have concluded that the 1/N expansion is 
most reliable on the BCS side of the crossover, let us 
now investigate this regime in more detail. The matrix 
of second derivatives of il'"' at the critical temperature 
needed to evaluate the NLO corrections /J^^^ and fi''^^ can 
be easily calculated from Eq. ([5]). One finds 



3nT 
4^' 



up to corrections which are exponentially small as kpa — > 
0— . (Note that we prefer to take the derivative with re- 
spect to A'^ instead of A since at the critical temperature 
we have to set A = afterwards.) The last needed coef- 
ficient, 



— (^-11 



is most easily obtained using the exact relation 
SirnT 



which follows from the observation that il'"^ can be ex- 
pressed in terms of a dimensionless function of the com- 
bination Pfi. 

The coefficient d^pfl^^'^ is strongly suppressed by the 
third power of the critical temperature so that it may be 
neglected and the required matrix of second derivatives 



[as in Eq. ([5]), just with modified variables] becomes 



1 



'A2/3 



-dA^M°^ 







As a result, the leading NLO contribution to the chemical 
potential decouples and is solely determined by Q^il'^' 



anda^^r!(o). 



(1) - 



This is not surprising since we know from Eq. (jlOp that 
in the BCS limit the chemical potential is governed by 
perturbative effects. The expression for the shift of the 
inverse critical temperature also simplifies to 



Pi''- 



dA^n^^^^ + 



dA^pnio) 



Substituting all the analytic expressions listed above as 
well as the chemical potential correction from Eq. pO|) 
we get the final prediction of the 1/N expansion for the 
NLO relative shift of the inverse temperature in the BCS 
limit, 



(1) 

C 

W 



4eF 
3n 



dA2n^'^^ + 



1 



9 -I- 2 log 2 
5^ 



kva 



where the only missing ingredient, that has to be evalu- 
ated numerically, is Oa^^^^^ 

The point of these considerations is that the slow, per- 
turbative convergence of the chemical potential results 
in a rather large offset in the critical temperature, which 
survives even in the limit kpa — > 0—. This is a con- 
sequence of the fact that c^As^il^''^ diverges in the BCS 
limit; we regard it an artifact of the expansion of the gap 
and number equations, leading to the expression ([4|. 

In addition to the critical temperature, we have also 
computed the gap and chemical potential at zero tem- 
perature, see Fig. [31 This has already been done by 
Veillette et al. although they did not investigate 

the asymptotic behavior in the BCS limit. Our calcula- 
tion differs from theirs only in that we have for technical 
reasons taken A^ instead of A as the variable to make 
the 1 /N expansion. Since the relative NLO correction to 
the gap at zero temperature is small, the effect of this 
change is nearly negligible. 

Fig. [3] suggests that the gap also acquires a constant 
offset that survives the BCS limit (even though we were 
not able to check this conclusion analytically). However, 
the offset is not the same as in the case of the critical tem- 
perature, as is clearly seen from Fig. |3]where we plot the 
ratio Tc/Aq. Thus, 1/N expansion predicts a departure 
of this ratio from the BCS value c'^/tt. In Sec. Illll we will 
present a calculation of the critical temperature within 
a 1/A^-inspired high-density approximation, where such 
artifacts will be absent. 
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FIG. 3: Gap and chemical potential at zero temperature as a 
function of the inverse scattering length. Dashed lines: MF 
approximation (LO in 1/N). Solid lines: NLO calculation 
based on expansion of A^. Dotted lines in lower panel: First- 
and second-order perturbative approximations to the chemi- 
cal potential, see Eq. (|10p . 



III. DENSE RELATIVISTIC MATTER 



—l/kpa 

FIG. 4: Ratio of the critical temperature to the gap at zero 
temperature as a function of the inverse scattering length. 
Dashed line: MF approximation. Solid line: NLO calculation 
with data taken from Figs. [1] and [S] Horizontal dotted line: 
Prediction of the BCS theory. 




FIG. 5: (a) RPA propagator of the collective mode in the case 
it transforms as an antisymmetric tensor of SU(A''), using the 
double-line notation, (b) One of the planar diagrams which 
dominate the thermodynamic potential in the large- A*' limit. 



We are now going to apply the 1 /N expansion to pair- 
ing in relativistic systems. Physically, this amounts to 
including the antiparticles among the degrees of freedom, 
and to modifying the fermion dispersion relation. In par- 
ticular the change of dispersion relation affects the ultra- 
violet structure of the theory, leading to new divergences 
that are not present in the nonrelativistic case [35|. In 
Sec. nil B I we evade this difficulty by making a suitable 
approximation, appropriate in the far BCS regime. 



A. NJL-type model 

Following closely the notation of Ref. [s^ , we consider 
a class of NJL-type models defined by the Lagrangian 

C = ^,{10 + /i7o - m)V' + I ^ H^lBQa^f, (11) 

a 

where ■0'' = C-ip'^ is the standard charge-conjugated 
Dirac spinor and the set of matrices Qa, acting on the in- 
ternal degrees of freedom, are normalized by Tr{QaQl,) = 
dab- Simplified as much as possible, this Lagrangian 



describes a system of interacting fermions with equal 
masses and chemical potentials. The pairing is assumed 
to occur in a spin-zero, positive-parity channel, but its 
flavor structure, determined by the matrices Qa, can be 
arbitrary. Once we understand in detail the fluctuation 
effects in this simple setting, we will move on to more 
realistic systems in our future work. 

The first task to solve is the generalization of the model 
to arbitrary so that we can subsequently make the ap- 
propriate expansion. With the application to quark mat- 
ter in mind, one may think it would be most convenient 
to use the color SU(3) symmetry already present in the 
system and extend it to SU(7V). We would like to explain 
here in detail why this would not work. 

As a consequence of the QCD interactions, quarks are 
assumed to pair in a color-antisymmetric configuration. 
(Symmetry would not make a difference. Important is 
that the Cooper pair cannot be a singlet of the sym- 
metry due to the complex nature of the group SU(3).) 
Upon the generalization of the theory the pairs would 
transform in the antisymmetric-tensor representation of 
SU(A^). Disregarding the fact that this would lead to a 
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very large number, of order 0{N'^), of collective modes, 
there is another, more serious problem. 

In Sec. nil we saw using the scaling of the thermody- 
namic potential that the 1/N expansion reproduces the 
MF approximation at the leading order. This may be 
also understood in terms of the collective mode propaga- 
tor. In the MF approximation, this consists of a geomet- 
ric series of graphs in the random-phase approximation 
(RPA). Now since the boson is a singlet of the Sp(2A^) 
symmetry, each fermion bubble contributes a factor N 
from the trace over the flavor space, thus compensating 
for factors 1/7V coming from the coupling g. The point is 
that when the boson becomes a tensor of SU(A^) rather 
than a singlet, the trace factors are lost, as is most eas- 
ily visualized with the help of the double-line notation, 
see Fig. [S](a). As a consequence, addition of each new 
fermion bubble in the RPA series suppresses the graph 
by another factor from the coupling so that the full 
series will not be resummed at any finite order in 

Is it then possible to make the coupling scale just as 
0(1)? No, because the thermodynamic potential is dom- 
inated by the planar diagrams with a single fermion loop 
(rather than the RPA ones), which contain the same 
number of flavor traces as the coupling factors, see Fig. [5] 
(b). So, in order that the power of A'^ in Feynman graphs 
is bounded from above and the 1/N expansion makes 
sense at all, the coupling has to decrease at least as 

We thus conclude that when the difermion field is de- 
fined to be a tensor of the SU(A^) color group, we will 
not get the RPA propagator and hence the pairing insta- 
bility at any finite order in The 1/A^ expansion in 
this case efficiently resums a different class of diagrams 
than necessary. To remedy this problem, we introduce 
a new quantum number to label the fermion fields, i.e., 
generalize Eq. to 

£ = + A*7o - m)i;, + j^Yl IVf TsOaV-^P- (12) 



This has two advantages. First, one does not need to rely 
on the presence of the color SU(3). The construction 
is absolutely general, regardless of the actual physical 
internal degrees of freedom. Second, the Cooper pair is 
a singlet with respect to the symmetry transformations 
acting on the new quantum number so that there are no 
additional bosonic degrees of freedom introduced by the 
extension of the symmetry and no unwanted NG bosons 
from its spontaneous breaking. On the other hand one 
has to set, just like in the nonrelativistic case, = 1 at 
the end of the calculation. 

With the above argument in mind, we proceed as in 
Sec. ini Upon bosonization of the theory p2)) by intro- 
ducing a set of auxiliary fields, 0a ~ ^i/JiJdQaipi, we 
arrive at the effective action 



S = n[ dr [ d^J-^^^^^^-NTT\ogg-'[M^,T)], 
Jo J 9 

(13) 

with the inverse fermion propagator in the Nambu space 
given by 



From the classical action we can generate the LO 
(RPA) propagator of the collective bosonic modes by a 
second functional derivative. Within this paper, we will 
for simplicity restrict our attention to the normal phase; 
the extension of the formalism below the critical temper- 
ature will be considered elsewhere. In the normal phase, 
the LO boson propagator becomes I?oafc = T^a^ab, 



1^-1/ X 1 1 

-Po (^^„,p) = - + - 



(2^ 



k , k 



/(ek^ + /i) + /(ek_ + m) - 1 , 1 - /(ek+ - m) - /(ek_ - /i)" 



2^ + ek+ + ek_ 



2/i - ek+ - ek_ 



1 - 



k 



Ck+Ck. 



/(ek+ + m) - /(ek_ - /i) , /(ek_ +Ai)-/(ek+ -/i) 



2/^ + Ck , - ek_ 



2^ + ek_ - Ck , 



(14) 



where k± = k ± ^ and f{x) = l/{e^^ + 1) is the Fermi- 
Dirac distribution function. 

Assuming that the order parameter fluctuations do not 
change the second order of the phase transition, we can 
find the critical temperature using the Thouless crite- 
rion [48j], i.e., by requiring that the normal-phase bo- 
son propagator has a pole (pairing singularity) at zero 
(four-)momentum. This is equivalent to the gap equation 
dn/dA^ = at A = 0. The expression dQ^^^dA'^ that 



appears in the NLO formula for the critical point (O, is 
thus seen to represent the one-loop boson self-energy at 
zero momentum. 

In a general scalar self-interacting theory the one-loop 
self-energy is given by the tadpole diagram with one quar- 
tic interaction vertex. In case of the theory defined by 
the action the effective four-boson vertex is gener- 
ated by a fermion loop, see Fig. [S) Let us concentrate 
on the flavor structure of the diagram. Assuming that 
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FIG. 6: One-boson-loop contribution to the collective mode 
self-energy. The double-dashed line denotes the resummed 
LO boson propagator. 



the external scalar legs carry the flavor indices a, h, the 
complete information about the pairing pattern will be 
encoded in the flavor trace 5cd TrF(QaQiQc(5d), the Kro- 
necker delta coming from the internal boson propagator. 
The combination 5cdQcQ\ ^^^^ quadratic Casimir oper- 
ator of the symmetry group in the representation of the 
Cooper pairs, and therefore must be proportional to the 
unit matrix as long as this representation is irreducible, 
i.e., we consider a single pairing channel [6l| . Writing 
5cdQcQ\ = C2I and taking the trace, we immediately 
find that C2 = ^c(i^c(i/dim((5), where dim((5) denotes 
the size of the matrices Qa- To conclude the argument 
we just note that this is equal to the number of internal 



fermionic degrees of freedom, iVp, while Scd^cd counts the 
number of bosonic degrees of freedom, Nb- The whole 
effect of the structure of the symmetry group will thus 
be the simple algebraic prefactor, 



ScdTTFiQaQlQcQi) 



Sab 



The full expression for the inverse boson propagator with 
one-loop correction then reads 

V-\P)^V,\P)+N^^dKdQVo{Q) 

X Tr [75^11 (/^ + P)l5Q22{K)-i5Qii{K + Qh5Q22{K)\ , 

where the subscripts n and 22 refer to the matrix struc- 
ture of the fermion propagator in the Nambu space. For 
the sake of brevity we used the notation for the four- 
momentum, P = (iujn,Tp), and the sum- integral. 



dK = T 



(2^ 



For zero external momentum the Matsubara sum in the 
fermion loop may easily be done and we arrive at the 
final analytic result, 



J 



D-i(0)=Po-i(0) 



/(a, 6; irtn 



DoHO) + N^pQVoiQ) E 
1 tanh^-Ktanh^-/3acosh-^^ ] 



1 + e/ 



m 



iQn + b + a 



1 tanh^ 



CkCk+q 

- tanh ^ 



•(k + q) 



a 



Iie^^J^i+^;^nn), (15) 



1 tanh ^ + tanh ^ 



with the usual relativistic notation, = ■\/k'^ + m^, 
= Ek + e/i. Note that upon taking the nonrelativistic 
limit, that is, including only particle degrees of freedom 
(with e = / = — ), shifting the chemical potential to 
Mnr = h — m, and approximating the dispersion rela- 
tion with e\t^ = m + ^ , we reproduce the expression for 
the one-loop propagator correction given by Nikolic and 
Sachdev 

B. High-density approximation 

The nonrelativistic limit of Eq. (|15p may be used di- 
rectly to calculate the NLO corrections to the critical 
temperature and chemical potential via Eq. ([5]) . To that 
end, we need to supplement Eq. (fTSj) with an analo- 
gous one-loop correction to particle density, which in turn 
yields the term (9^17'"'^' in Eq. This is how the results 
presented in Fig. [T]were obtained. 

On the contrary, in the full relativistic description the 



one-loop term in Eq. is badly divergent. This is, 

of course, no surprise since in a relativistic scalar self- 
interacting field theory the one-loop tadpole graph has a 
quadratic divergence. Nevertheless, this divergence has 
nothing in common with the many-body physics, and can 
be removed by renormalizing the parameters of the the- 
ory in the vacuum. In order to avoid this complication 
and also the interference of all energy scales from the 
high-energy vacuum physics down to the scale of Cooper 
pairing, we resort to a high-density approximation [49|, 
which is appropriate in the far BCS region where the 
pairing energy scale is well separated from the Fermi 
scale. This approximation is at the MF level known to 
soften the ultraviolet divergences and give some cutoff- 
independent predictions such as the universal BCS ratio 
of the gap at zero temperature and the critical tempera- 
ture [50|. 

In our one-loop calculation we have to be more careful. 
We therefore spell explicitly all simplifying assumptions 
that we make. First, we neglect antiparticle contribu- 
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tions. This can be appropriate even in the ultrarelativis- 
tic Hmit as long as the pairing gap/critical temperature is 
much smaher than the Fermi energy so that the relevant 
excitations are the quasiparticles and quasiholes near the 
Fermi surface. In a strongly coupled relativistic super- 
conductor, the antiparticle effects may be non-negligible 
po| but will not change our conclusions qualitatively. 

Second, we approximate, as usual, the volume measure 
for the integral over fermionic momentum by 

where ^ = is the energy with respect to the Fermi level 
and J\f = ixk-p/^Tz^ is the density of states on the Fermi 
surface. This brings in one subtlety. In a nonrenormaliz- 
able theory such as the NJL model, there is an inherent 
ambiguity in the way we label internal propagators of the 
Feynman graphs with momenta which satisfy momentum 
conservation in the interaction vertices. For instance, in 
the one-loop pairing susceptibility one often labels 
the propagators with k, k ± p instead of k + ^ , k — ■! 
used here, and imposes a cutoff on the loop momentum 
k. This ambiguity can be removed by a proper renor- 
malization which makes the graph finite. However, once 
we make the above introduced replacement of the inte- 
gration measure in the high-density approximation, the 
new measure is no longer translationally invariant. That 
is, the result depends on the momentum assignment to 
the propagators even in an otherwise finite loop integral. 
In Eq. (|14p we chose the symmetric momentum assign- 
ment because it reduces cutoff dependence of the result 
and also leads to a physically intuitive suppression of the 
high-momentum pair modes by Pauli blocking. 

Third and most importantly, in order to retain in the 
calculation only the physically relevant degrees of free- 
dom, we cut off the ^ integration at the pairing scale, 
1^1 < A. In quark matter, such a cutoff is effectively intro- 
duced in terms of the momentum-dependent gap function 
by solving the QCD Schwinger-Dyson equations [Slj . In 
addition, a cutoff much smaller than the Fermi energy is 
necessary in order to make the high-density approxima- 
tion consistent |4§|. In practice, we verify that the results 
are not sensitive to a precise value of the cutoff by doing 
all calculations for two different values, A = 2ri°-' and 

Fourth, once we restrict ourselves to the low-energy ex- 
citations about the Fermi surface, we expand the disper- 
sion relations in terms of the Fermi velocity and an "effec- 
tive mass" . This allows us to introduce efficiently dimen- 
sionless variables and treat on the same footing the non- 
relativistic (NR) limit of Sec. [ITl as well as the opposite- 
extreme, ultrarelativistic (UR) limit of zero fcrmion mass 
(which is a reasonable approximation for quark matter 
composed solely of the u and d flavors). Concretely, tak- 
ing the LO critical temperature TjP^ as a unit for the 
energy variables and ri^^/wp, where is the Fermi ve- 
locity, as a unit for external momentum, the dispersions 



in the two opposite limits acquire very similar forms, 
^k+p = <^k + IpI cose* H- — NR limit, 

4 ^ 

Ck+p = Ck + IpI cos 6* H — sin^ 9, UR hmit, 

2 

where ^ = ^/Ti"^ and p — pwf/Tc"', is the angle be- 
tween the vectors k and p, and the nonrelativistic chemi- 
cal potential is understood in the first line, without spec- 
ifying the subscript nr in the following. Using this pro- 
cedure, all integrals would factorize into a product of 
powers of the temperature and Fermi velocity and a uni- 
versal dimensionless function of the ratio TjP^ / were it 
not for the coupling g. 

The last step in the construction therefore has to be 
the renormalization of the bare coupling. In the context 
of atomic gases near unitarity it is customary to do this 
by fixing the s-wave scattering length at zero momentum 
in the vacuum as in Eq. However, in an effective 

description near the Fermi surface, this is no longer con- 
venient. We therefore renormalize the bare coupling with 
the help of the gap equation at zero temperature, which 
has the same divergence structure as the inverse propa- 
gator HH). This is effectively done by the replacement 

i-/(ekj-/(ekj 

it^n - ^k+ - Ck_ 

, i-/(aj-/(ek_) ^ 1 

^ ii^n - a+ - Ck_ 2y/(^ + A2 

in the inverse propagator (|14p in the high-density ap- 
proximation. Now everything is expressed in terms of 
the dimensionless ratios Tc^^ /fi and Ti^V^o- The ratio 
Tc"' / ^J■ is used as the input parameter which measures the 
strength of the interaction. The ratio Tc^^ / Aq is equal 
to the BCS value — « 0.567 in the infinite-cutoff limit. 

7r 

With the explicit cutoff on the ^-integration, we adjust 
the value of Ag appropriately in order to ensure that the 
Goldstone theorem is satisfied and the propagator (flU) 
has an exactly massless pole. 

Due to the explicit cutoff, the loop part of the in- 
verse propagator (jl4p drops rapidly for external momenta 
larger than the cutoff as a result of Pauli blocking. The 
boson propagator approaches a constant value, equal to 
g/N. This is natural: At large momentum, pairing fluc- 
tuations are suppressed and the auxiliary field propagator 
recovers the original contact four-fermion interaction. In 
order that we really include just the effect of the fluctua- 
tions of the order parameter, we make in the second term 
of Eq. (fT5|) the replacement 

In terms of Feynman graphs, this means removing from 
the RPA series of diagrams the first, constant term, keep- 
ing all other terms that involve multiple rescattering of 
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FIG. 7: Sum of the LO and NLO contributions to the inverse 
boson propagator. The asymptotically constant part of the 
fermion self-energy in the second diagram may be absorbed 
in a perturbative renormalization of the Fermi energy. 




FIG. 8: Correction to the critical temperature calculated 
using the high-density approximation in the nonrelattvistic 
limit. The calculation was done with the cutoff on fermion 
energy variable set to 2ri''' (solid line) and AT^:^^ (dashed 
line). 



the two fermions in the pair. Formally, this subtraction 
can be justified by observing that the diagram in Fig. 
ini may also be viewed as a fermion loop with the inser- 
tion of a one-loop fermion self-energy, see Fig. [71 The 
fermion propagator with the insertion of g /N is nothing 
else than the first-order perturbative correction by the 
contact four-fermion interaction. We can then pick this 
constant term out of the one-boson-loop diagram and 
add it to the one-fermion-loop MF graph contributing to 
Dq^(0), where it is absorbed in the perturbative renor- 
malization of the Fermi energy |62| . 



C. Results 

Recalling finally from the BCS theory that at temper- 
ature T different from the MF critical temperature ri^-* , 
Vq^{0) is equal to AAlog(r/Ti°^), setting ^'-^(O) to zero 
leads to the following compact formula for the critical 
temperature, 



exp 



Nf 



(0) 




10 



10-1 , , 10-2 

y(0) 



10 



-3 



(16) 



FIG. 9: Correction to the critical temperature calculated 
using the high-density approximation in the ultrarelativistic 
limit. The calculation was done with the cutoff on fermion 
energy variable set to 2ri''' (solid line) and 4ri'^' (dashed 
line) . 



where X is a dimensionless function (different in the two 
limits), given by the integral in Eq. p5p in the dimen- 
sionless variables introduced above. 

Eq. pB)) constitutes our final result regarding relativis- 
tic superconductors, written in the most general form: 
The coupling constant is completely eliminated in favor 
of the MF ratio TjP'^ / ^. Also, the specific form of the 
pairing channel only enters through the algebraic factor 
Nb/Nf. 

Two comments are in order here. First, Eq. (|16p gives 
the correction to the critical temperature at fixed chemi- 
cal potential; we do not solve the number equation along 
with the Thouless criterion to obtain a result at fixed 
density. In the context of color-superconducting quark 
matter, the (baryon number) chemical potential is usu- 
ally treated as a free parameter. In fact, it is a more 
suitable parameter than the density itself in the case that 
the phase diagram involves first-order phase transitions, 
where the density becomes discontinuous. 

Second, the 1/iV algorithm based on Eq. ([5]) (reduced 
to a one-variable problem by fixing the chemical poten- 
tial) would suggest to interpret the exponent in Eq. (fTB]) 
as the relative change of the critical temperature, or mi- 
nus the relative change of the inverse critical tempera- 
ture, depending on the choice of variable. In this section, 
we used the 1/A^ expansion to derive Eq. (jlSp as the 
one-loop corrected Thouless criterion, and to evaluate 
the loop correction at the MF value of the critical tem- 
perature. We now go slightly beyond the 1/N philosophy 
in the sense that our result, Eq. p^ . does not need any 
further expansion and defines the critical temperature Tc 
in an unambiguous manner. The value of thus cal- 
culated is always positive, no matter how large the loop 
function X is. (Of course, we would still consider the 
used approximation unreliable once X becomes of order 
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one or larger.) In addition, the temperature correction 
does not display a finite offset in the BCS limit, as found 
in Sec. ITTl and assigned to the 1/A^ expansion as an ar- 
tifact. On the contrary, it drops rapidly with decreasing 
ratio Tc^^ as one would naively expect [63j . 

To complete the discussion of the results, we show in 
Figs. [5] and [5] the numerically calculated values of the 
function X in the NR and UR limits, respectively. For 
Tc^^ / ^ smaller than about 0.1 the functions may be very 
well approximated by a simple empirical power law, 

(17) 

which can be used for a fast rough estimate of the size of 
fluctuation effects. 



D. Possible impact on QCD phase diagram 

Finally, we wish to illustrate a possible impact of fluc- 
tuations on the QCD phase diagram. In Fig. [10] we 
display the phase diagrams from a simple NJL-model cal- 
culation; the model is the same one as adopted in Ref. 
[S^ l with the diquark coupling Gd chosen such that the 
CFL gap is 50MeV at = 500 MeV in the chiral SU(3) 
limit. Fig. [TO] (a) shows the phase diagram at vanishing 
strange quark mass, i.e., in the chiral SU(3) limit. The 
dashed line is the critical temperature at the leading or- 
der in 1/A'' expansion, and the suppression of the critical 
temperature due to fluctuation effects at NLO is included 
by means of the analytic formula p6p . Interestingly, be- 
cause N-q/N-p = 1 in the CFL case while it is just 1/2 in 
the 2SC phase, we expect a finite region with 2SC pairing 
below the normal phase even in the chiral SU(3) limit. 

However, since the approximation which led to Eq. 
(I16p is only valid in the high-density regime 

where ri^V^ 

is small, we have to keep in mind that the estimates are 
not quantitative at low density. It is also important to 
note that we derived the shift of critical temperature of 
the CFL phase, taking into account the fluctuations in 
the normal phase. Now that we know that the 2SC phase 
interposes between the normal and CFL phases, it would 
be more appropriate to somehow take into account the 
fluctuations within the 2SC phase for a more realistic es- 
timate of the temperature of the phase transition between 
the CFL and 2SC phases. 

Fig. [TO](b) shows the phase diagram at a flnite strange 
quark mass, but without quantum fluctuation effects. 
The strange quark mass is set to Ms — 200 MeV, and 
is for simplicity treated as an external parameter rather 
than a dynamical one. Also, the intricate charge neu- 
trality effects are ignored here because they are not im- 
portant for our purposes; we just remark that the charge 
neutrality constraints bring fine splittings to the T^s re- 
sulting in the appearance of tiny regions with dSC or 
uSC pairing [H, [s^. Comparing the two plots, we can 



see that the quantum fluctuation effects may play a sig- 
nificant role at high density which may be similar to that 
of the strange quark mass. This can be understood in 
a model- independent way as follows. Within a weak- 
coupling Ginzburg-Landau approach, it was shown that 
the strange quark mass and charge neutrality result in 
shifts in the melting temperatures Ti of the order param- 
eters Ai of the CFL phase of the order [s^ 

i^^-Miogf^^ 

These corrections die rapidly as (M^/^)^ at large ^, while 
the correction to T^"'' from order parameter fluctuations 

prevails as it behaves asymptotically like Tc^^ / fj,. Con- 
sidering in addition the fact that Ms is a decreasing func- 
tion of ^ whereas the superconducting gap turns out to 
increase as ^ oo [H, ]55[ , we conclude that the quan- 
tum fluctuation is more important than the effect of a 
strange quark mass at high density. 

In a more realistic situation, there would be an- 
other source of fluctuations from thermal configurations 
of gauge fields [S^ which makes the superconducting- 
normal transition first order. The shift of the critical 
temperature turns out to be positive and is estimated 
to be proportional to the QCD coupling g in the weak- 
coupling regime [56|. Therefore, in realistic quark mat- 
ter the order parameter and gauge field fluctuations will 
compete each other. This issue certainly deserves a fur- 
ther study in future. 

IV. SUMMARY AND CONCLUSIONS 

We have investigated the 1/N expansion for strongly 
interacting Fermi systems, proposed recently [ll,[l3|- We 
flrst studied in detail the case of nonrelativistic Fermi 
gas near unitarity, extending the previous results by the 
calculation of the critical temperature off the unitarity. 
Even though the 1/A^ expansion can give a result for 
the critical temperature at unitarity which is very close 
to the prediction from Monte Carlo simulations 
there is an inherent ambiguity due to the choice of ob- 
servable to generate (and truncate) the 1/A'' series. This 
ambiguity makes the 1/A'' expansion in the current set- 
ting useless when the corrections to the MF theory are 
large, in particular on the BEC side of the crossover. 

We paid particular attention to the evolution of the 
fluctuation corrections in the BCS regime, where they 
are expected to be small, and the 1/A^ series thus to con- 
verge fast. We showed that the next-to-leading order in 
the 1 /N expansion reproduces the well-known perturba- 
tive correction to the chemical potential up to second 
order. As far as the critical temperature is concerned, 
the fluctuation correction indeed decreases at weak cou- 
pling, but leaves a flnite offset in the kpa — > 0— limit. 
We argued that this is likely to be an artifact of the 1 /N 
expansion. 
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FIG. 10: (a) Phase diagram at vanishing strange quark mass Ms = but with quantum fluctuation efi^ects. (b) The same phase 
diagram for Ms — 200 MeV without fluctuation efl^ects. ri"'' is the critical temperature at the leading order in 1/A^. 



In Sec. mil we applied the idea to strongly coupled rel- 
ativistic superconductors, having in mind in particular 
color-superconducting dense quark matter. We used a 
simple class of NJL-type models and resorted to a high- 
density approximation in order to avoid conceptual diffi- 
culties associated with renormalization and entanglement 
of several energy scales. Our results are summarized in 
Eqs. (HH) and ((T7I) and Figures [H] and O They are physi- 
cally intuitive in the sense that the fluctuation corrections 
are small in the BCS limit, decreasing linearly with the 
ratio Tc^^ / fi, i.e., exponentially with the inverse coupling, 

and become large as xi'^^ / ^ approaches the order of 0.1 
and further grows. 

In particular for typical color superconductors, the cor- 
rections to critical temperature are expected to be as 
large as tens percent. Another important conclusion is 
that, within the simple setting used here, the fluctuation 
corrections are expressed as a universal function of the 
dimensionless ratio Tc^"^ /^. The whole dependence on 
the symmetry structure of the pairing is encoded in an 
algebraic prefactor, which counts the number of bosonic 
and fermionic degrees of freedom. This makes it straight- 
forward to compare the effects for different competing 
superconducting phases, and thus estimate the impact of 
the order parameter fluctuations on the phase diagram 
such as in Fig. [TOl Upon this investigation of the fluctu- 



ation corrections within a simple model, we plan, in our 
future work, to include the effects of chemical potential 
mismatch and color neutrality in order to obtain a more 
realistic description of quark matter. 

Finally, we would like to stress the conceptual simplic- 
ity of this approach to order parameter fluctuations. We 
do not need to solve a complicated set of self-consistent 
integral equations like in other techniques going beyond 
the MF approximation, such as the Cornwall-Jackiw- 
Tomboulis one. Instead, one just has to evaluate a single 
multidimensional sum-integral. We therefore believe that 
the 1/A^ expansion may provide an efficient tool to de- 
termine the fluctuation effects in such strongly-coupled 
systems as, for instance, the color superconductors. 
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It was argued in Ref. [iq ] that due to the distinct physical 
origin of the LO (fermions) and NLO (bosons) contribu- 
tion to the thermodynamic potential, their relative size 
should not be used to test the accuracy of the 1/N ex- 
pansion. This postpones the issue of convergence of the 
series to higher orders, yet the ambiguity encountered at 
NLO remains. 

Strictly speaking, our calculation differs from the orig- 
inal NSR theory [3| in that we perform the Matsubara 
sum in the one-loop contribution to the number density 
directly while it is more conventional to make analytic 
continuation to the real frequencies and use spectral rep- 
resentation. Of course, the results should coincide as long 
as the analytic continuation has been done properly. 
Should we include several pairing channels, our results 
would naturally generalize, giving a different fluctuation 
effect for each of them even in the normal phase. 
Actually, the first-order perturbative renormalization of 
the fermion propagators in the fermion loop in T>q^{0) 
produces two equal terms, one from each propagator. 
However, this is exactly what we need because the loop 
diagram, being symmetric, carries a symmetry factor i. 
As shown by Gorkov and Melik-Barkhudarov, particle- 
hole fluctuations lead to a suppression of t he g ap and 
critical temperature even in the BCS limit [57l.l58| . We do 
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not expect this effect to be reproduced by our cafcufation 
since it inciudes only pair ffuctuations [45l |. 



